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Abstract 

Background: Expectation maximizing (EM) is one of the common approaches for image segmentation. 

Methods: an improvement of the EM algorithm is proposed and its effectiveness for MRI brain image 
segmentation is investigated. In order to improve EM performance, the proposed algorithms incorporates 
neighbourhood information into the clustering process. At first, average image is obtained as neighbourhood 
information and then it is incorporated in clustering process. Also, as an option, user-interaction is used to improve 
segmentation results. Simulated and real MR volumes are used to compare the efficiency of the proposed 
improvement with the existing neighbourhood based extension for EM and FCM. 

Results: the findings show that the proposed algorithm produces higher similarity index. 

Conclusions: experiments demonstrate the effectiveness of the proposed algorithm in compare to other existing 
algorithms on various noise levels. 
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1. Background 

The application of image processing techniques for medical 
imaging process rapidly increases. Most medical images are 
stored and represented in softcopy [1]. Ultrasound, X-ray 
computed tomography, digital mammography and mag- 
netic resonance imaging (MRI) are the most common 
medical imaging types [2], MRI can give different grey level 
for different tissues and various types of neuropathology if 
its acquisition parameters are adjusted [3]. 

Data acquisition, processing and visualization techni- 
ques facilitate diagnosis. Medical image segmentation 
plays a very important role in many computer-aided 
diagnostic tools. These tools could save clinicians' time 
by simplifying complex time-consuming processes [4]. 
The main part of these tools is to design an efficient seg- 
mentation algorithm. Medical images mostly contain 
unknown noise [5], in- homogeneity [6] and complicated 
structures. Therefore, segmentation of medical images is 
a challenging and complex task. Medical image segmen- 
tation has been an active research area for a long time. 
There are many segmentation algorithms but there is no 
generic algorithm for a totally successful segmentation of 
medical images [7] . 
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Clustering methods are common for MRI brain segmen- 
tation. Expectation-maximization (EM) and fuzzy c-mean 
(FCM) are the most popular clustering algorithms. The 
Gaussian mixture model (GMM) is a popular segmenta- 
tion method. EM is used to estimate the parameters of 
this model. FCM and EM only consider the intensity of 
images and in noisy images, intensity is not trustful [8-10]. 
Usually, spatially adjacent pixels belong to the same clus- 
ter. Many algorithms introduced to make FCM [11-17] 
and EM robust against noise but nevertheless most of 
them were and are flawless to some extent. Usually, spa- 
tially adjacent pixels belong to the same cluster. Many 
researchers attempted to incorporate spatial information 
into FCM and EM to overcome the noise problem. Zhang 
et. al. [18] proposed a novel Gaussian hidden Markov 
Random Field (HMRF) model to integrate spatial informa- 
tion into Gaussian model. They used a Markov Random 
Field-Maximum A Posteriori (MRF-MAP) approach to 
estimate the model solution. Recently, Tang et al. [19] 
proposed a neighbourhood-weighted Gaussian mixture 
model to overcome misclassification on the boundaries 
and on inhomogeneous regions of MRI brain images with 
noise. A. R. F. d. Silva [20] proposed two Bayesian algo- 
rithms (DPM, rjMCMC) which use Markov chain sam- 
pling techniques to find normal mixture models with an 
unknown number of components. They used algorithms 
for MRI segmentation and compared performance of their 
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algorithms with published results for two exist Bayesian 
based MRI brain segmentation methods (KVL [21], MPM- 
MAP [22]). 

Gonzalez Ballester et al. [23] and Tohka et al. [24] 
reported a statistical models namely a novel trimmed 
minimum covariance determinant (TMCD) for the esti- 
mation of the parameters of partial volume models to 
address partial volume averaging. 

In order to make Gaussian mixture model more robust 
against complex tissue spatial layout, Greenspan et al. [25] 
proposed the parameter-tied, constrained Gaussian mix- 
ture model (CGMM) to capture this problem. The mix- 
ture model composed of a large number of Gaussians for 
each tissue is used to capture the complex tissue spatial 
layout. The Gaussian parameters of a tissue are tied using 
intensity as global feature. The parameters are learned 
using the expectation- maximization (EM) algorithm. 

In [26], a nonparametric Bayesian model, known as 
Dirichlet process mixture model (DPMM) is proposed to 
overcome the limitations of current parametric finite mix- 
ture models. The DPMM permits unknown number of 
components in the mixture and allow robust segmentation 
of brain with unknown or incomplete specifications. 

In [27], local cooperative unified segmentation (LOCUS) 
approach based on distributed local MRF models for brain 
segmentation is presented. The volume is partitioned into 
sub volumes and a set of local and cooperative Markov 
random field (MRF) models are distributed. In order to 
ensure consistency, neighbour local MRFs are estimated 
cooperatively. The intensity in-homogeneity correction is 
not required due to precisely fit of Local estimation with 
the local intensity distribution. 

In this paper, a new modification to GMM and EM is 
introduced by incorporating neighbourhood information 
into likelihood function and EM steps. The average of 
neighbour pixels around each pixel is calculated prior to 
GMM clustering and incorporated in GMM and EM func- 
tions beside the pixel value. 

The rest of this paper is organized as follows. The stan- 
dard GMM model and EM segmentation algorithm are 
presented in Section 2.1. In Section 2.2, proposed modified 
EM algorithm is described. Also, improvement of segmen- 
tation results using use-interaction is presented in section 
2.3. Experimental and comparison results are presented in 
Section 3 and this paper is concluded in Section 4. 

2. Methods 

A modification to GMM is introduced by incorporating 
neighbourhood information into likelihood function and 
EM steps. 



(voxel) with M mixing coefficients. Each component is 
assigned to one target class and the goal is to obtain the 
class probabilities of each pixel (voxel). The probability 
distribution of the yth component is denoted by Pj(Xi\9j), 
where x t is pixel i in input image and 6j is the parameter 
(mean ftj and covariance matrix £/) of component /. The 
probability distribution each pixel (voxel) can be 
described as a mixture of probability distributions as 
follows: 



p{xi\6) = J2 a iPj( x i\°i) 
1 



(1) 



ydet(27r£ ; ) 



Where a ; denotes the mixture coefficient with the 

M 

constraint, J2 a i = 1 The probability distribution of 

i=i 

component ; is modelled by a Gaussian distribution 
with mean ftj and covariance matrixZ ; : 



pj{Xi | Oj) =pj(Xi | fAj,Zj) 



(2) 



Usually, maximum likelihood (ML) estimation is used 
to find the parameters. The log-likelihood expression for 
the parameter 6 and the image X is defined as follows: 



log(L(0 |X))=logf]K*i|0) 



i=i 



N M 



(3) 



Finding the ML solution from this equation is difficult. 
Usually, the expectation-maximization (EM) is used to 
obtain the parameters. EM steps are demonstrated in 
the following: 

E-step. Bayes' rule is used to obtain the probability of 
data Xi belong to class 6j (E-step): 



®iPi[xi I Of) 



M 



E«/M*»l*/) 



(4) 



M-step. Probability obtained in E-step is used to 
obtain mixing coefficient, mean and covariance matrix 
(M-step): 



2.1 Standard GMM 

The Gaussian mixture model assumes M mixed compo- 
nent densities (Gaussian distribution) for each pixel 



1 N 

< 1 = ^E^u^ £ ) 
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c. EM steps are repeated until convergence. 
2.2. Modified GMM 

The average of neighbour pixels around X{ is calculated 
prior to GMM clustering. In the likelihood function 
(Equation 3), distribution value of X{ is added to the dis- 
tribution value of pixel Xi as neighbourhood information: 



\o%{L{0\X))=\ogY[p{x l \0) = 



i=i 



(8) 



J>g£>|[(l - P) * Pj {xi | Of) + p*Pj{xi | Of)]) 
i=i j=i 

The parameter /3 determines the weight of neighbour- 
hood information. Incorporating neighbourhood infor- 
mation improves the performance of segmentation 
methods in high level of noise, but the blurring effect 
degrades the performance of them in low noise level. In 
order to overcome the degrading effect of algorithms in 
low level of noise, the variance of noise is used to spe- 
cify the weight of neighbourhood information Its 
value is set to O", where a is the variance of noise. In 
previous neighbourhood based EM extensions, neigh- 
bourhood information is calculated in clustering itera- 
tion; but in this algorithm X{ is computed before 
iteration, thus, the clustering will be faster. An extension 
of EM named EM-1 is introduced to solve likelihood 
function. The EM is modified as follows: 

a. In Equation 4, distribution value of X{ is added to the 
distribution value of pixel x t as neighbourhood informa- 
tion: 

A=[{l-p)*p j {x i \0f) + p*p j {xi\0j)\ 

a r A (9) 



p{j\x i/ e t ) 



E4 A 



b. In Equation 6, X{ is added to x t as neighbourhood 
information: 



/z 



J2{{l-P)*x l + p*x i )p(j\x u 6 t ) 

t+1 - - (10) 



c. In Equation 7, the distance of X{ from the compo- 
nent centre is added to the distance of x t from the com- 
ponent centre as neighbourhood information: 



d(x) = (x-f^ +l )(x-^ +l ) 



J2p(j\x i ,0 t Ud{x i )+p.d{x i )) 



t+i _ i=i 



(11) 



1=1 



In MRI, noise behaves as Rician distributed noise. 
Rician noise approaches Gaussian distribution in high 
Signal to Noise Ratio (SNR) and Rayleigh distribution in 
low SNR [28]. Rician distribution in the background is 
Rayleigh because there is no signal. The Rayleigh PDF 
of the statistically independent observations is 



P({0,}) = [1 7T2 e 



-(Of)/(2a 2 



(12) 



i=i 



Where O is observations and cr^ oise is the variance of 
noise. The variance of noise is obtained by maximizing 
the log-likelihood of PDF function with respect to var- 
iance: 



Noise 



1 H 

i=l 



(13) 



In other words, background pixels are considered as 
observations (O) and the variance of noise is obtained 
applying equation 13 on background pixels values. For 
that, the powers of background pixels values are com- 
puted and half of the average of resulted values is con- 
sidered as variance of the noise. 

Also, in-homogeneity correction [6] is applied to input 
image with in-homogeneity pollution and the propose 
GMM is applied on in-homogeneity corrected image. 

2.3. Improving Segmentation Results Using User 
Interaction 

Sometimes, due to in-homogeneity, low contrast, noise 
and inequality of content with semantic, automatic meth- 
ods fail to segment image correctly. Therefore, for these 
images, it is necessary to use user interaction to correct 
method's error [29]. However, robust semi-automatic 
methods can be developed in which user interaction is 
minimized. 
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Sometimes, segmented image, for example in Figure 1 
(b), either has pixels from two or more tissues in one clus- 
ter (csf and grey matter of brain in cluster number 2) or 
pixels from one tissue in two or more clusters (white mat- 
ter in clusters number 2 and 3). For solving this problem, 
user selects clusters contain several tissues (cluster num- 
ber 1) to be re-clustered to two sub clusters. Figure 1(c) 
demonstrates sub clusters of class number 1. The cluster 
number 1 is clustered to two sub clusters number 11 
and 12. 

This process continues until user is satisfied. That 
means quality of segmentation depends on user. Then, 
to solve problem of several clusters for one tissue, user 
selects clusters for each tissue (clusters 12 is also 
selected for grey matter). Steps of this method listed as 
follow: 

1. Input volume is clustered to the n clusters where 
n is the number of target class (tissues). The output 
is clustered volume. 

2. Under segmentation: If some clusters contain more 
than one target class (tissue), user selects such clus- 
ters to be partitioned more; each user selected cluster 
is re-clustered to two sub clusters. This process con- 
tinues till user is satisfied. The output is clustered 
volume without under segmentation. 

3. Over segmentation: If several clusters correspond to 
one target class (tissue), user selects clusters for each 
target class. The output is final clustered volume. 



3. Experimental Results and Discussion 

The proposed extension of EM (EM-1) and the existing 
neighbourhood-based extension of EM [19] (referred as 
NWEM in this paper for clear understanding) are simu- 
lated and tested on the simulated volumes from BrainWeb 
[30] and real volumes from Internet Brain Segmentation 
Repository (IBSR) [31], 



Moreover, reported results on simulated volumes for 
existing extensions of EM (DPM, rjMCMC, KVL, MPM- 
MAP) and existing neighbourhood based extension for 
FCM (FCM_S [32], FCM_EN [33], FGFCM [34], FLICM 
[35] and NonlocalFCM [36]) are used to evaluate pro- 
posed algorithm. 

Also, the reported results on real volumes from IBSR 
are used to evaluate proposed algorithms. Furthermore, 
mentioned FCM extensions simulated and tested on real 
volumes. 

The results of algorithms are compared quantitatively 
to analyse their performance. The neighbourhood size, N 
for proposed algorithm is set to 3 x 3. Three indices 
(similarity index, false positive ratio and false negative 
ratio) [37] are used to evaluate the algorithms quantita- 
tively. The similarity index p t of class i is the degree of 
the class pixels matching between ground truth and seg- 
mentation result for the same class. The false positive 
ratio rfp represents extra pixels of class i and the false 
negative ratio rfn represents lost pixels of class i. They 
are defined as follows: 



2 I Xi n Yi | 

+ \ Yi I 



rfpi 



\Yj\-\ Xi n Yj I 
\Xt\ 



rfrii = 



I K | - | Xi n Yj | 
I ^ | 



(14) 



where X t represents class i in ground truth and Y t 
represents the same class in the segmentation result. 
Each index for full segmentation results is the average 
of that index for all classes. 

3.1. Simulated volumes 

The simulated MRI volumes are obtained from Brain- 
Web. A simulated data volume with Tl-weighted 
sequence, slice thickness of 1 mm and a volume size of 
217 x 181 x 181 is used. Non-brain tissues are removed 
prior to segmentation. 

The number of tissue classes in the segmentation is 
set to three: grey matter (GM), white matter (WM) and 
cerebrospinal fluid (CSF). All pixels in the image are 
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contributed in segmentation process but in evaluation 
process, background pixels are ignored following pre- 
vious works utilized in this paper. In the public data- 
bases which have been used in the paper and generally 
in brain MRI volumes, background pixels have black 
value. Therefore, cluster with lowest average grey value 
is considered as background. 

First, EM-1 and NWEM were applied to a slice of Tl- 
weighted brain image corrupted by different noise levels. 
Figure 2 and Figure 3 show the segmentation results of 
applying the afore-mentioned algorithms on a Tl-weighted 
normal brain slice in the presence of 9% and 7% rician 
noise, respectively. 

The segmentation results of white matter (WM), grey 
matter (GM) and cerebrospinal fluid (CSF) are depicted 
in. (a) is noisy image, (b) is ground-truth, (c) to (d) are the 
segmentation results of NWEM and EMI, respectively. 

From the above qualitative comparison, it was not dif- 
ficult to find that NWEM was more influenced by the 



noise in comparison with EMI, in which fewer artefacts 
were evident, resulting in clearer segmentation result. 

Also, the proposed segmentation algorithm (EM-1) 
and NWEM are applied to brain volume and average 
similarity value is used to evaluate them. Figure 4 shows 
the average similarity indexes p of mentioned algorithms 
in different noise levels. 

Figure 4 shows that EM-1 produces higher similarity 
indexes and lower rfp and rfn, meaning that this algo- 
rithm produces more accurate segmentation results. The 
similarity index of EM-1 decreases more slowly than 
NWEM algorithm when noise level increases. In the 
same time, the rfp and rfn of EM-1 increases faster than 
NWEM algorithm. 

Both algorithms give similar results, under 5% noise 
level. However, for more than 5%, EM-1 exhibits much 
better results than the NWEM algorithm. Incorporating 
average of neighbourhood information, in clustering 
process of NWEM, make this algorithm robust against 




a) 



c) 



b) 




d) 



Figure 2 The segmentation results of applying EMI and NWEM on a slice of image with 9% Rician noise, (a) Noisy image, (b) Ground- 
truth, Segmentation results of (c) NWEM and (d) EM1. 
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Figure 3 The segmentation results of applying EMI and NWEM on a slice of image with 7% Rician noise, (a) Noisy image, (b) Ground- 
truth, Segmentation results of (c) NWEM and (d) EM1. 



noise but has blurring as side effect. It seems that with 
increasing noise level more than 5% noise level; this 
incorporation cannot overcome high level of noise. 

Also the effect of different neighbourhood sizes on 
performance of proposed segmentation algorithm (EM- 
1) is investigated. Figure 5 shows the average similarity 
index p of EM-1 for different neighbourhood sizes on 
volume with 9% noise. Figure 5 shows that when the 
neighbourhood size is increased, the similarity index of 
EM-1 decreases sharply. This means blurring effect in 
EM-1 depends on neighbourhood size. 

The speed of EMI and NWEM in segmenting a slice 
was also investigated. Figure 6 represents the average 
time required to segment a slice using the mentioned 
algorithms. Figure 6 shows that EMI is faster than 
NWEM. The neighbourhood information in NWEM is 
calculated in NWEM clustering iteration. Therefore, it is 
time-consuming. 

The proposed segmentation algorithm is also compared 
with current extensions for EM. The average similarity 



indexes p for proposed algorithm (EM-1) and several 
current extensions for EM (DPM, rjMCMC, KVL and 
MPM-MAP) are shown in Figure 7. Figure 7 shows that 
EM-1 produces highest similarity indexes. The proposed 
segmentation algorithm gives results comparable with 
the best reported results, in low level of noise. However, 
for noise levels more than 5%, EM-1 algorithm outper- 
form other competing algorithms and this difference in 
performance gets more in 9% noise level. 

Also, EM-1 is compared with current existing neigh- 
bourhood based extensions for FCM. Figure 8 shows the 
average similarity indexes p for EM-1 and FCM extensions 
(FCM_S, FCM_S1, FCM_EN, FGFCM and FLICM) in dif- 
ferent noise levels. At 3% noise level, the results for pro- 
posed segmentation algorithm and the best reported result 
were close. Above 3% noise, EM-1 produces higher simi- 
larity index and were the most convincing in segmenta- 
tion. The superiority of these algorithms increases with 
increasing in noise level. FLICM shows worst performance 
it seems it is not suit algorithm for brain volumes. 
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Figure 4 The average similarity indices p, rfp and rfn for different noise levels. 
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In [25], the parameter-tied, constrained Gaussian mix- 
ture model (CGMM) is applied on image volume from 
brainweb with different noise levels. Average similarity 
index for different algorithms with variant noise levels 
(3%, 5%, 7%, 9%) are: CGMM (0.93, 0.93, 0.92 and 0.895) 
and KVL (0.925, 0.915, 0.895 and 0.865). The proposed 
segmentation algorithm outperforms KVL and CGMM. 

3.2. Real volumes 

The superiority of our algorithm is also demonstrated on 
real MRI volumes. The real MRI volumes are obtained 
from the IBSR by the Centre for Morphometric Analysis 
at Massachusetts General Hospital. 20 normal data 



volume with Tl-weighted sequence are used. First, pro- 
posed algorithm (EM-1) is applied to slices of a real MRI 
volume with size 256*256*53. The average similarity 
index p for volume image is 0.7986. Figure 9 shows the 
similarity indexes of proposed algorithm (EM-1) for each 
slice of MRI volume. In almost all slices, the proposed 
algorithms exhibit better results for WM in compare to 
results for GM. Better performance of proposed algo- 
rithms in WM is due to more simplicity and compactness 
of WM in compare to GM. 

EM-1 and NWEM are applied to all 20 normal real MRI 
volumes and average similarity index p is used to compare 
the segmentation results, quantitatively. Figure 10 shows 
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Figure 5 The average similarity index p for different neighbourhood sizes on simulated volume with 9% noise. 



the average similarity index, rfp and rfn values of both 
algorithms for all 20 normal volumes. Figure 10 shows 
that EM-1 outperforms NWEM. EM-1 produces higher 
average similarity indexes p and lower rfp and rfn. 

The average similarity index values of proposed algo- 
rithm for 20 normal real MRI volumes and EM exten- 
sions (reported results in IBSR) are compared. Figure 11 
shows the average similarity index values of different 
algorithms for all 20 normal volumes. Figure 11 shows 
that the similarity index for proposed segmentation 
algorithms is higher than competing methods. It can be 
seen clearly that proposed algorithm has a better perfor- 
mance over reported results, meaning that proposed 
algorithm produces more accurate segmentation results 

The proposed algorithms are also compared with 
neighbourhood based extensions for FCM. Figure 12 
shows the average similarity indexes p for proposed 
algorithm and FCM extensions (FCM_S1, FCM_EN, 
FGFCM) for all 20 normal volumes. 

It can be seen clearly that proposed algorithm has a 
better performance over FCM extension methods, and 



produces more accurate segmentation results. FCM 
extensions also incorporate neighbourhood information 
in FCM clustering process, but, it seems that incorporat- 
ing neighbourhood information improves EM more than 
FCM method. 

In [24], a novel trimmed minimum covariance deter- 
minant (TMCD) method an extension for Gaussian mix- 
ture model is applied on 20 normal image volumes from 
IBSR. The average jaccard value was 0.6722. The average 
jaccard values for EM-1 is: 0.695. The similarity index 
for EM-1 is higher than reported result, meaning that 
EM-1 produces more accurate segmentation results. 

In [25], the parameter-tied, constrained Gaussian 
mixture model (CGMM) is applied on 18 volumes 
from 20 normal image volumes (except volume 4-8 
and 202-3) in IBSR website. The CGMM results is 
compared with reported results from the IBSR website, 
as well as with the Marroquin algorithm [38]. Marro- 
quin's algorithm is an atlas-based Bayesian segmenta- 
tion algorithm. The CGMM algorithm outperforms 
other studied methods. Jacc similarity index CGMM 




Algorithms 



Figure 6 Average times required to segment a slice using the proposed algorithm (EMI) and NWEM. 
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Figure 7 The average similarity indices p for different noise levels. 
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Figure 8 The average similarity indices p for EM-1 and FCM extensions in different noise level. 
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Figure 9 The similarity index of proposed algorithm when applied for real volume. 






Figure 10 The average similarity index, rfp and rfn of proposed algorithm when applied on 20 real volumes. 
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Figure 1 1 The average similarity index of different algorithms when applied on 20 real volumes. 




was: 0.67. The average jaccard values for EM-1 is: 
0.6971. EM-1 outperforms the best reported result 
which is for CGMM. 

In [26], a nonparametric Bayesian model, known as 
Dirichlet process mixture model (DPMM) is applied on 13 
volumes (1_24, 2_4> 5_8, 6_10, 7_8, 11_3, 12_3, 13_3, 
15_3, 16_3, 100_23, 110_3 112_2) from the 20 normal Tl- 
weighted brain image volumes from IBSR. The similarity 
index for DPMM is higher than competing methods. Dice 
similarity index for DPMM was: 0.7071. The proposed 
algorithms are applied on the same volumes. The average 
Die value for EM-1 is: 0.8219. The similarity index for pro- 
posed method is higher than the best reported result 
which is for DPMM, meaning that proposed method are 
the most convincing in segmentation. 



In [27], local cooperative unified segmentation 
(LOCUS) approach which is based on distributed local 
MRF models for brain segmentation is applied on the 
20 normal Tl -weighted brain image volumes from IBSR. 
LOCUS-T is compared with published results for SPM5 
and FAST. Die similarity index for different methods 
are: LOCUS-T = 0.765, SPM5 = 0.81, FAST = 0.765. 
The average Die value for EM-1 is: 0.8211. EM-1 out- 
performs the best reported result which is for SPM5. 

Also, improvement of segmentation result using user- 
interaction is investigated. Proposed algorithms and the 
same algorithm with user-interaction are applied to all 20 
normal real MRI volumes and similarity index p is used 
to compare the segmentation results, quantitatively. The 
average similarity index values of both algorithms for 
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Figure 12 The similarity index of proposed algorithm and neighbourhood based FCM extensions when applied on 20 real volumes. 



different volume were presented in Figure 13 and Table 
1. Figure 13 shows that user-interaction improves perfor- 
mance of proposed algorithm and increases similarity 
indexes p in all image volumes. 

4. Conclusion 

In this paper, an extension of EM has been introduced. In 
order to overcome the problem of standard EM in the pre- 
sence of noise, the introduced algorithms are formulated 
by modifying the equations of the standard EM algorithm 
which allow the neighbourhood pixels to be incorporated 
in the labelling of a pixel. Introduced algorithm is tested 



on simulated MRI volumes, with different noise levels and 
real volumes. The performance of the existing neighbour- 
hood based EM and FCM algorithms and proposed algo- 
rithm are compared qualitatively. 

The similarity index, p is used to evaluate different 
algorithms. Experiments demonstrate the effectiveness 
of the proposed algorithm in compare to other existing 
algorithms on various noise levels in terms of similarity 
index, p. 

In future, we consider doing research on other kinds 
of segmentation methods to improve their functional- 
ities. Also, we will analyse the effects of different 
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Figure 13 The similarity index of different algorithms when applied on 20 real volumes. 
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Table 1 The similarity index of different algorithms when applied on 20 real volumes 



Volumes Algorithms 


EM-1 with user-interaction 
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0.8377 



clustering methods in segmentation of medical images 
for the diagnosis of abnormal or various important mat- 
ters in medical images. 
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